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Using a simple analytic formalism, we demonstrate that significant dark matter self-interactions 
produce halo cores that obey scaling relations nearly independent of the underlying particle physics 
parameters such as the annihilation cross section and the mass of the dark matter particle. For dwarf 
galaxies, we predict that the core density p c and the core radius r c should obey p c r c fts 41 Mopc“ 2 
with a weak mass dependence ~ M 0,2 . Remarkably, such a scaling relation has recently been 
empirically inferred. Scaling relations involving core mass, core radius, and core velocity dispersion 
are predicted and agree well with observational data. By calibrating against numerical simulations, 
we predict the scatter in these relations and find them to be in excellent agreement with existing 
data. Future observations can test our predictions for different halo masses and redshifts. 


Introduction. The standard cosmological model with 
a cosmological constant A and cold dark matter (CDM) 
has been confirmed on large scales by a wealth of suc¬ 
cessful predictions. Yet despite its simplicity and exper¬ 
imental success, a variety of challenges to the ACDM 
paradigm persist. On the observational side, the CDM 
paradigm suffers from the well-known cusp-core problem 
(see |Tj for a review); numerical simulations HU! predict 
that the density of dark matter halos should rise sharply 
towards the center of the halo forming a “cusp”, whereas 
observations (see e.g. EH7] and references therein) in¬ 
dicate that the dark matter in the center of halos rises 
more gently, resembling a “core” even in dwarf galax¬ 
ies where the baryonic content is negligible [8]. Moti¬ 
vated by observations, various models of interacting dark 
matter have been proposed to alleviate this tension in¬ 
cluding scattering HUD] and annihilation HU, although 
some authors, e.g [T2], argue that dark matter cores do 
not appear upon a more careful analysis. Many inter¬ 
esting experiments (for recent progress see e.g. [33J E] 
and references therein) have been devised to constrain 
some combination of mass and interaction cross section. 
Here, we focus on the model-independent observational 
signatures of annihilation. Remarkably, we find that scal¬ 
ing relations exist that are generic predictions of self¬ 
interactions where any dependence on the cross section a 
and the particle mass cancel. We can thus make predic¬ 
tions about annihilation signature without introducing 
any additional parameters beyond what is required in a 
minimal ACDM cosmology. Throughout this Letter, we 
assume h = ^/(lOOkms^ 1 Mpc -1 ) = 0.7, f l m = 0.27, 
and n A = 0.73 US]. 

Core scaling relations To derive the core scaling re¬ 
lations, we need two simple ingredients. The first ingre¬ 
dient is the initial condition of the dark matter density 
profile, before enough cosmic time has passed for annihi¬ 
lation to significantly alter the profile. We take the dark 
matter profiles to be of the Navarro-Frenk-White (NFW) 
form H: 


p(r,t = t 0 ) 


_ Po _ 

{r/r s ){l + r/r s ) 2 ' 


( 1 ) 


Here to is the age of the universe when the dark matter 
halo virialized. We will comment shortly on the modifi¬ 
cation of our results if the inner profile of the dark matter 
halo is more accurately described by the Einasto profile 
[IB] , which is likely a better fit to simulations [TTJ . To 
determine po and r s for a halo of mass M and redshift 
z, we adopt an expression for the concentration param¬ 
eter C 2 oo(M,z) = r s /?' 200 j where ?' 2 oo is defined as the 
radius interior to which the average density is 200 times 
the critical density of the universe p C rit = (3H 2 /8nG), 
where H is the Hubble parameter. 

We adopt the fitting formulae of Dutton and Maccio 
118], their equations 7, 10, and 11, calibrated from the 
results of their N-body simulations: 

log c 2 oo = a + 61og(M/10 12 /i _1 M Q ) 
a = 0.520 + 0.385 exp (-0.617Z 1 - 21 ) (2) 

b = —0.101 + 0.026.3. 


These formulae are in agreement with a recently proposed 
universal model of PH- 

The second ingredient is the evolution of the dark mat¬ 
ter profiles with time. We will only consider the simplest 
case of s-wave annihilation where (av) = const. This 
is a reasonable approximation for most annihilation pro¬ 
cesses unless a quantum selection rule prevents s-wave 
annihilation or the cross section acquires a sharp veloc¬ 
ity dependence due to Sommerfeld enhancement. Ignor¬ 
ing the gravitational back reaction, the time derivative 
of the density is given by 

p(x,t) =-Tp 2 (x,t), (3) 


where Tp has units of inverse time. The solution to this 
equation is 


P(r, t) 


Pir, t 0 ) 

1 + fTtp{r, t o y 


( 4 ) 


where we have inserted a fudge factor / > 1 (effec¬ 
tively rescaling T) in anticipation that gravitational back- 
reaction should change the profile. (If gravitational back- 
reaction is negligible, / = 1). 
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Plugging in our initial condition yields 


P(r,t) 


_ Pc _ 

(r/r c )( 1 +r/r s ) 2 + 1’ 


( 5 ) 


where we define a core density and core radius p c = 
(/T t)^ 1 and r c = (/ T t p 0 ) r s , respectively. For the mo¬ 
ment, we simply define r c the core radius, we will show 
shortly that it corresponds approximately to the notion 
of core radius defined by other authors. 

Since different authors define these quantities differ¬ 
ently, it is worth noting what these quantities represent, 
so that they can be sensibly compared to observational 
work. It is straightforward to verify that both quantities 
can be expressed in terms of properties of the central den¬ 
sity profile: the central density is simply p c = p(r = 0), 
and Taylor expanding the logarithmic slope 


dlogp 

dlogr 



shows that r c is the length scale over which the logarith¬ 
mic slope significantly deviates from 0. A full compari¬ 
son of the behavior of our density profiles with existing 
work is given in Figure 1, which shows that the numerical 
value of r c (obtained by fitting to data) may differ from 
the values of (r c ) Bu rkert or (r c ) Za vaia, depending on the 
precise value of r s . However, the difference is unlikely to 
be greater than ~ 0.1 to 0.2 log increments, which trans¬ 
lates to a difference of several tens of percent. In other 
words, for a large range in possible values of r a , our def¬ 
inition of r c is consistent with (r c ) Bu rkert or (r c ) Z avai a up 
to a factor of order unity. 

These non-parametric properties of the central den¬ 
sity profile allow us to connect our definitions with other 
parametric profiles. For example, the above considera¬ 
tions lead us to identify (p 0 ) Burkert = Pc and (r c ) ?urkert « 
r c where the subscripts denote definitions found in Burk- 
ert 0, since the core radii in both our profile and Burk- 
ert’s profile are turnover radii, e.g., the length scale over 
which the logarithmic slope becomes non-negligible. 

Note that we have not made any assumptions about 
the constancy of / with particle physics or halo parame¬ 
ters. Kaplingliat et al. Eli] have argued that 3 < / < 10 
(decreasing with mass) accurately accounts for the grav¬ 
itational back-reaction in the limit that annihilation oc¬ 
curs slowly such that the adiabatic invariant rM(r) is 
conserved; numerical simulations of annihilating HUES 
and self-interacting (scattering only) [2Tl dark matter 
confirm the accuracy of this approximation. These simu¬ 
lations also show that our results are qualitatively unaf¬ 
fected by the choice of NFW or Einasto profiles, at least 
when the dark matter interaction is strong enough so that 
the core radius r c > 10 -2 r s and the differences between 
NFW and Einasto are significant. 

More quantitatively, Zavala et al. [25j showed that self¬ 
interacting dark matter profiles with a/m = 0.1 cm 2 g^ 1 
are well fit by a profile of the form (using our notation) 


PZavala(^) 


_ Pc _ 

(r/r c )( 1 + r/r s ) 2 + (1 + r/r s ) 2 ' 


( 7 ) 


For small radii r <C r s , the last term (l + r/r s ) 2 ss 1 
and we recover equation ([5]). For large radii r r s , 
the r 3 term in the denominator will dominate, and the 
quadratic term (1 + r/r s ) 2 can be neglected. In either 
regime, our model is in agreement with the results of 
Zavala et al. [22:. The typical volume averaged difference 
between our models (p — p Za vaia)/po is not more than a 
few percent within 2 r s , and is entirely negligible for much 
larger volumes. 

Although we have motivated equation (3) by consid¬ 
ering annihilating dark matter, we note that scattering 
also has a similar effect, as scattering will kick parti¬ 
cles from the high density inner regions to larger radii, 
where their contribution to the mass budget is negligible. 
Since the scattering rate is also proportional to p 2 , T in 
equation (3) should also include contributions from the 
scattering cross section. These results are in agreement 
with numerical simulations (e.g. the results of Elbert 
et al. [51]) which predict a flat slope interior to a core 
radius ~ r c and a profile which converges to the orig¬ 
inal NFW form at large radii. Note, however, that if 
a self-interacting dark matter halo is in the gravother- 
mal collapse regime, our results do not apply. Indeed 
numerical simulations [53J [53] show that during collapse, 
the inner density profile becomes significantly more cuspy 
than a constant density core. The failure of our formal¬ 
ism in this regime is not surprising, as processes such 
as gravothermal collapse obviously do not preserve adia¬ 
batic invariants. 

Although r c and p c will have a large spread due to the 
different formation times of halos, there exists an overall 
scaling relation 


p c r c = Por s . ( 8 ) 

We emphasize that the above equation was not derived by 
making any assumptions about whether or not r c <C r s . 
Note that the formation time t , the fudge factor /, and 
the annihilation rate T no longer appear in this result for 
the core surface density. (However if the cross section is 
negligible, r c —> 0 and the above scaling relation will be 
impossible to observe). The quantities on the right hand 
side of equation (6) are properties of the dark matter halo 
which can be easily computed. It is a straightforward to 
integrate equation (1) to find 

= jy t j,„( fe i,) 2/3 MV> 

log(C200 + 1 ) — C200 /(1 + C200) 

where the normalization N = 10(10/7r) 1 / 3 /(3 2 / 3 ) ss 7.07. 
Since C 2 qo(M 1 z) is a very weak function of mass, we ex¬ 
pect that the right hand side to scale roughly with M 1 / 3 . 
A more detailed calculation shows that the mass depen¬ 
dence is even weaker por s oc M 0 2 at low redshifts and 
even at a relatively high redshift z = 3, por s oc M °- 3 . 
Thus, if we consider a subclass of halos such as dwarf 
galaxies, p c r c will be approximately constant. This 
is precisely the scaling relation inferred empirically in 
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log (r/r c ) 

FIG. 1: Comparison of our density profile with Burkert [7] and Zavala et al. [22] , On the y-axis is the logarithmic 
slope of the density profile, dp/d\ogr , on the x-axis is the logarithmic radius in units of r c . The gray lines show our 
profile for r s = 4 r c to r s = 16 r c in log increments. Similarly, the pink dashed lines show different values of r s using 
the parameterization found in Zavala et. al. The thick black line is from [7]. Note that observations provide the best 
constraints in the region r/r c < 1. This graph illustrates that our definition of r c is approximately consistent with 
other definitions of r c found in the literature. For example, the value of r c inferred by fitting a density profile to 
some given data may be somewhat larger or smaller compared to the value of (r c )Burkert or (r c )z a vaia- 




FIG. 2: Predicted p c r c as a function of mass M 200 and redshift 2 . As sample points, we display the Milky Way [25] 
in black, a typical dwarf galaxy in gray, and the Phoenix galaxy cluster [2~6j in red. On the right panel, different 
curves display different halo masses M 200 , which increase in unit log-increments. On the left panel, the different 
curves represent different redshifts which increase in increments of 0.5. 


Spano et al. [J, Donato et al. [5], Kormendy and Free¬ 
man [ 6 ] and most recently Burkert [7j. 

Ultimately, the near-constancy of p$r c is a consequence 
of the fact that the initial conditions provide by (non¬ 
interacting) LCDM simulations have a density profile 
that scales approximately oc r _1 . (The resulting profiles 
with interaction taken into account do not need p oc r^ 1 ). 
For example, if dark matter halos were more accurately 
described by older secondary infall models 123, the inner 
profiles would be isothermal: p oc r -2 , which would not 
lead to the near-constancy of p c r c . 

Comparison with observations. Although it is obser- 
vationally difficult to measure the corresponding M 2 00 of 
a dwarf galaxy, numerical simulations [281 suggest that a 
typical value is M 200 ~ 10 1 o Mq. Using this as a fiducial 
value, we predict 


p c r c = 41 Mq pc 2 x 


M 2 00 
1O 1O M 0 


0.18 


( 10 ) 


where the logarithmic slope ( M b with b = 0.18) has been 
obtained by appropriately (log-log) linearizing the mass 
dependence at the fiducial value. Although the slope is 
calibrated at the fiducial value, it remains approximately 
constant over a large range of masses; even for a galaxy 
cluster b = 0.22. In particular, for the Phoenix galaxy 
cluster (25], p c r c w 1.1 x 10 3 M 0 /pc -2 . We show the full 
dependence on mass and redshift in Figure 2. 

Remarkably, such a scaling relation has been empiri¬ 
cally inferred by Burkert )7], who reports that (in our no¬ 
tation) p c r c = 75+45 A/© pc -2 over 18 magnitudes in blue 
magnitude, covering a sample that ranges from dwarf 
galaxies to giant galaxies. Here we have identified the 
central density p(r = 0 ) = p c and the core radius r c as 
the length scale associated with a turnover in logarithmic 
slope. The reported uncertainties are not la uncertain¬ 
ties but encompass all but 1 or 2 of the 48 data points. Al¬ 
though the median dark matter halo mass of their sample 
is not easily measured and therefore not reported, a rea¬ 
sonable value is M = 10 11 ' 5 Mq which lies between dwarf 
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galaxies (M ~ 10 1o Mq) and giant galaxies (M ~ 10 13 ). 
Taking M = 1O 115 M 0 and 2 = 0, we have from equation 
(7) that p c r c = 78A/ 0 pc -2 , which fully consistent with 
the value that is empirically inferred. 

The predicted scaling of p c r c with mass can also be 
tested. Since the luminosity L* of a galaxy is propor¬ 
tional to the number (and thus the mass Af* of stars, 
L* oc M*. For low mass galaxies M 20 o < 1O 121 M 0 , ob¬ 
servations [29] indicate M 2 oo cx M* with /? = 0.33lg 3 5 , 
with the slope gradually steepening to oc M 9 ' 75 for very 
massive galaxies. Hence, for low mass M 2 oo < 1O 121 M 0 
(intermediate mass M 200 ~ 1O 121 M 0 ) galaxies, we pre¬ 
dict that p c r c oc M 2 5q 8 oc M 306 oc Lj with 7 = 0.06 
(7 = 0 . 10 ), values which are consistent with the observed 
slope in [ 6 |: p c r c oc Lj, where 7 = 0.058 ± 0.067. Future 
observations, particularly of a sample of massive galaxies 
where the M* - L* relation is steeper, could further con¬ 
strain 7 in order to provide a more stringent test of our 
predictions. We also note that our model predicts virtu¬ 
ally no redshift evolution for 2 < 1. Hence a detection of 
redshift dependence in p c r c would falsify our model. 

It is worth mentioning that our derivation of the scal¬ 
ing relations allows us to compute the expected scatter 
in the observed relations. The two sources of scatter in 
our model are the scatter in c 2 oo that exists even for 
a fixed mass and redshift, and the scatter in M if the 
observed population of dwarf galaxies contains galaxies 
of different virial masses. Note that from equation (7), 
p c r c oc c 200 for large values of c 2 oo, whereas the explicit 
mass dependence is much weaker. A typical amount of 
scatter (at the ler level) associated with c 2 oo is ~ 0.1 
to ~ 0.2 dex [18]. Unless the standard deviation in the 
virial masses of an observed population of dwarf galaxies 
is greater than ~ 1 to 2 orders of magnitude, the scatter 
in c 2 oo will dominate. To compare to empirical results, 
we will assume that the population of dwarf spheroidal 
galaxies observed in 0 does not contain such a diversity 
of masses. Certainly if the scatter in mass were much 
larger than an order of magnitude, the fact that all of 
their halos have log Mo .3 = 7.1 ± 0.3 where Mo .3 is the 
mass enclosed in the inner 0.3 kpc would seem peculiar, 
since the total mass is expected to be a strong function 
of Mo. 3 . In particular, for an NFW profile gives leads 
to a mass dependence [3'D] M 20 o oc (M 0 . 3 ) 2 ' 9 , though de¬ 
viations from p oc r _1 in the inner regions will change 
this dependency. Assuming that c 2 oo has a ler scatter 
of 0.15 dex, we predict using Monte Carlo methods that 
p c r c = 78+23 -^0 pc -2 with a 68% confidence interval. 
Considering that the uncertainties reported in Burkert 
[7: contain ~ 96% of the data and are thus close to 2cr 
bounds, there is good agreement between the predicted 
and observed scatter. 

Our scaling laws are also in agreement with the older 
results of Spano et al. @j and Donato et al. [5], which 


find in their sample of galaxies and dwarf galaxies p c r c = 
lO 2 'i 5±o ’ 2 M 0 pc -2 [5]. The somewhat higher value (by a 
factor of ~ 2) of p c r c compared to the results of Burkert 
[7] could be due to the fact that their sample is dominated 
by more massive halos. 

As a consistency check, we may derive additional scal¬ 
ing relations from our model. For example, if we assume 
that r c <C r s (an assumption that we so far have not yet 
made use of), the density in the region r < r c is approx¬ 
imately constant, and 

Me 4 tt 

— ~ -ttPcTc ~ const, (11) 

r c o 

and the core velocity dispersion er c obeys ~ 

GM(r < r c )/r 2 ~ Gpor c . Hence, 

cr 2 r“ 1 ~ const. (12) 

The last two scaling laws have also been empirically in¬ 
ferred in [7] and have been shown to be consistent with 
available data. However, we caution that these additional 
results can only be derived in our formalism if r < r e , 
which is not the general case. 

Conclusion. Starting from an NFW profile and a sim¬ 
ple treatment of annihilation, we have derived a universal 
scaling relation p c r c ~ 41 M 0 pc -2 x (Af 2o o/1O 1O M 0 ) 0 ' 18 
at 2 = 0. Said differently, it has been remarked that 
anomalies with ACDM are associated with an accelera¬ 
tion scale ao ~ 10 -9 cms -2 [5] where new physics be¬ 
comes relevant. In units where c = G = 1, a character¬ 
istic acceleration is equivalent to a characteristic surface 
density. In this work, we have derived the surface density 
scale (see equation 9), assuming that the new physics is 
dark matter self-interactions. 

Remarkably, the derivation did not involve any addi¬ 
tional parameters beyond what is needed in LCDM. This 
scaling relation holds independently of the dark matter 
mass and annihilation or scattering cross sections as well 
as the amount of adiabatic expansion / experienced by 
the core. Our results are also relatively insensitive to 
slow motions of baryonic matter in the dark matter halo, 
provided that the effect of the baryons is to adiabatically 
expand or contract the core. This scaling relation is thus 
a robust signature of self-interaction, independent of the 
detailed properties of the dark matter particles. Special¬ 
izing to the case of dwarf galaxies, we have shown that 
both the magnitude of p c r c and the scatter are in ex¬ 
cellent agreement with the data. We have also checked 
the model by deriving additional scaling relations which 
also agree with experiment. The predicted evolution of 
halo properties with redshift and mass can be tested with 
future observations. 
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